AI delivers Michaelis constants as fuel for genome-scale metabolic models

Michaelis constants (Km) are essential to predict the catalytic rate of enzymes, but are not widely available. A new study in PLOS Biology uses artificial intelligence (AI) to accurately predict Km on a proteome-wide scale, paving the way for dynamic, genome-wide modeling of metabolism

Multistationary and Oscillatory Modes of Free Radicals Generation by the Mitochondrial Respiratory Chain Revealed by a Bifurcation Analysis

The mitochondrial electron transport chain transforms energy satisfying cellular demand and generates reactive oxygen species (ROS) that act as metabolic signals or destructive factors. Therefore, knowledge of the possible modes and bifurcations of electron transport that affect ROS signaling provides insight into the interrelationship of mitochondrial respiration with cellular metabolism. Here, a bifurcation analysis of a sequence of the electron transport chain models of increasing complexity was used to analyze the contribution of individual components to the modes of respiratory chain behavior. Our algorithm constructed models as large systems of ordinary differential equations describing the time evolution of the distribution of redox states of the respiratory complexes. The most complete model of the respiratory chain and linked metabolic reactions predicted that condensed mitochondria produce more ROS at low succinate concentration and less ROS at high succinate levels than swelled mitochondria. This prediction was validated by measuring ROS production under various swelling conditions. A numerical bifurcation analysis revealed qualitatively different types of multistationary behavior and sustained oscillations in the parameter space near a region that was previously found to describe the behavior of isolated mitochondria. The oscillations in transmembrane potential and ROS generation, observed in living cells were reproduced in the model that includes interaction of respiratory complexes with the reactions of TCA cycle. Whereas multistationarity is an internal characteristic of the respiratory chain, the functional link of respiration with central metabolism creates oscillations, which can be understood as a means of auto-regulation of cell metabolism. © 2012 Selivanov et al

Unveiling a key role of oxaloacetate-glutamate interaction in regulation of respiration and ROS generation in nonsynaptic brain mitochondria using a kinetic model.

Glutamate plays diverse roles in neuronal cells, affecting cell energetics and reactive oxygen species (ROS) generation. These roles are especially vital for neuronal cells, which deal with high amounts of glutamate as a neurotransmitter. Our analysis explored neuronal glutamate implication in cellular energy metabolism and ROS generation, using a kinetic model that simulates electron transport details in respiratory complexes, linked ROS generation and metabolic reactions. The analysis focused on the fact that glutamate attenuates complex II inhibition by oxaloacetate, stimulating the latter's transformation into aspartate. Such a mechanism of complex II activation by glutamate could cause almost complete reduction of ubiquinone and deficiency of oxidized form (Q), which closes the main stream of electron transport and opens a way to massive ROS generating transfer in complex III from semiquinone radicals to molecular oxygen. In this way, under low workload, glutamate triggers the respiratory chain (RC) into a different steady state characterized by high ROS generation rate. The observed stepwise dependence of ROS generation on glutamate concentration experimentally validated this prediction. However, glutamate's attenuation of oxaloacetate's inhibition accelerates electron transport under high workload. Glutamate-oxaloacetate interaction in complex II regulation underlies the observed effects of uncouplers and inhibitors and acceleration of Ca2+ uptake. Thus, this theoretical analysis uncovered the previously unknown roles of oxaloacetate as a regulator of ROS generation and glutamate as a modifier of this regulation. The model predicted that this mechanism of complex II activation by glutamate might be operative in situ and responsible for excitotoxicity. Spatial-time gradients of synthesized hydrogen peroxide concentration, calculated in the reaction-diffusion model with convection under a non-uniform local approximation of nervous tissue, have shown that overproduction of H2O2 in a cell causes excess of its level in neighbor cells

Integration of enzyme kinetic models and isotopomer distribution analysis for studies of in situ cell operation

A current trend in neuroscience research is the use of stable isotope tracers in order to address metabolic processes in vivo. The tracers produce a huge number of metabolite forms that differ according to the number and position of labeled isotopes in the carbon skeleton (isotopomers) and such a large variety makes the analysis of isotopomer data highly complex. On the other hand, this multiplicity of forms does provide sufficient information to address cell operation in vivo. By the end of last millennium, a number of tools have been developed for estimation of metabolic flux profile from any possible isotopomer distribution data. However, although well elaborated, these tools were limited to steady state analysis, and the obtained set of fluxes remained disconnected from their biochemical context. In this review we focus on a new numerical analytical approach that integrates kinetic and metabolic flux analysis. The related computational algorithm estimates the dynamic flux based on the time-dependent distribution of all possible isotopomers of metabolic pathway intermediates that are generated from a labeled substrate. The new algorithm connects specific tracer data with enzyme kinetic characteristics, thereby extending the amount of data available for analysis: it uses enzyme kinetic data to estimate the flux profile, and vice versa, for the kinetic analysis it uses in vivo tracer data to reveal the biochemical basis of the estimated metabolic fluxes

3D UWB magnitude-combined tomographic imaging for biomedical applications. Algorithm validation

Biomedical microwave imaging is a topic of continuous research for its potential in different areas especially in breast cancer detection. In this paper, 3D UWB Magnitude-Combined tomographic algorithm is assessed for this recurrent application, but also for a more challenging one such as brain stroke detection. With the UWB Magnitude-Combined concept, the algorithm can take advantage of both the efficiency of Fourier Diffraction Theorem-based tomographic formulation and the robustness and image quality improvement provided by a multi-frequency combination.Peer ReviewedPostprint (published version

UWB High-contrast robust tomographic imaging for medical applications

In this paper a complete UWB Circular Tomographic System robust to high contrast or large objects, applied to Breast Tumor Detection, is presented. The main contribution of this paper is to focus on the implementation of a two degrees of freedom imaging setup in order to deal with non-symmetric objects and to demonstrate its functionality with realistic breast phantoms.Peer ReviewedPostprint (published version

Edelfosine-induced metabolic changes in cancer cells that precede the overproduction of reactive oxygen species and apoptosis

Background: Metabolic flux profiling based on the analysis of distribution of stable isotope tracer in metabolites is an important method widely used in cancer research to understand the regulation of cell metabolism and elaborate new therapeutic strategies. Recently, we developed software Isodyn, which extends the methodology of kinetic modeling to the analysis of isotopic isomer distribution for the evaluation of cellular metabolic flux profile under relevant conditions. This tool can be applied to reveal the metabolic effect of proapoptotic drug edelfosine in leukemia Jurkat cell line, uncovering the mechanisms of induction of apoptosis in cancer cells. Results: The study of 13C distribution of Jukat cells exposed to low edelfosine concentration, which induces apoptosis in ¿5% of cells, revealed metabolic changes previous to the development of apoptotic program. Specifically, it was found that low dose of edelfosine stimulates the TCA cycle. These metabolic perturbations were coupled with an increase of nucleic acid synthesis de novo, which indicates acceleration of biosynthetic and reparative processes. The further increase of the TCA cycle fluxes, when higher doses of drug applied, eventually enhance reactive oxygen species (ROS) production and trigger apoptotic program. Conclusion: The application of Isodyn to the analysis of mechanism of edelfosine-induced apoptosis revealed primary drug-induced metabolic changes, which are important for the subsequent initiation of apoptotic program. Initiation of such metabolic changes could be exploited in anticancer therapy

Targeted Metabolomics Shows That the Level of Glutamine, Kynurenine, Acyl-Carnitines and Lysophosphatidylcholines Is Significantly Increased in the Aqueous Humor of Glaucoma Patients

The composition of the aqueous humor of patients with glaucoma is relevant to understand the underlying causes of the pathology. Information on the concentration of metabolites and small molecules in the aqueous humor of healthy subjects is limited. Among the causes of the limitations is the lack of healthy controls since, until recently, they were not surgically intervened; therefore, the aqueous humor of patients operated for cataract was used as a reference. Sixteen aqueous humor samples from healthy subjects undergoing refractive surgery and eight samples from glaucoma patients were used to assess the concentration of 188 compounds using chromatography and mass spectrometry. The concentration of 80 of the 188 was found to be reliable, allowing comparison of data from the two groups (glaucoma and control). The pattern found in the controls is similar to, but not the same as, that reported using samples from “controls” undergoing cataract surgery. Comparing data from glaucoma patients and healthy subjects, 57 of the 80 compounds were significantly (p < 0.05) altered in the aqueous humor. Kynurenine and glutamine, but not glutamate, were significantly increased in the glaucoma samples. Furthermore, 10 compounds were selected considering a statistical score of p < 0.0001 and the degree of change of more than double or less than half. The level of C10 (decanoyl)-carnitine decreased, while the concentration of spermidine and various acyl-carnitines and lysophosphatidylcholines increased in glaucoma. Principal component analysis showed complete segregation of controls and cases using the data for the 10 selected compounds. The receiver operating characteristic curve these 10 compounds and for glutamine allowed finding cut-off values and significant sensitivity and specificity scores. The concentration of small metabolites in the aqueous humor of glaucoma patients is altered even when they take medication and are well controlled. The imbalance affects membrane components, especially those of the mitochondria, suggesting that mitochondrial abnormalities are a cause or consequence of glaucoma. The increase in glutamine in glaucoma is also relevant because it could be a means of keeping the concentration of glutamate under control, thus avoiding its potential to induce the death of neurons and retinal cells. Equally notable was the increase in kynurenine, which is essential in the metabolism of nicotine adenine dinucleotides

Toward 3D UWB tomographie imaging system for breast tumor detection

A novel 3D tomographic algorithm for short range cylindrical geometries using UWB frequency range is presented. The algorithm has been applied to breast tumor detection, nevertheless, its non application-specific character permits the use in other applications. The detection capability of the tomographic algorithm is proved through numerical simulations and experimental measurements of canonical and more realistic bodyattached breast phantoms. For the acquisition of the experimental data, a virtual array-based cylindrical measurement setup has been built. At the current stage of development 3 degrees of freedom are available allowing the scanning of bodies with revolution symmetry.Peer ReviewedPostprint (published version